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0^ , A stochastic variational inequality is proposed to model an elasto-plastic oscillator excited by a filtered 

white noise. We prove the ergodic properties of the process and characterize the corresponding invari- 
ant measure. This extends Bensoussan-Turi's method (Degenerate Dirichlet Problems Related to the 
Invariant Measure of Elasto-Plastic Oscillators, AMO, 2008) with a significant additional difficulty of 
' increasing the dimension. Two points boundary value problem in dimension 1 is replaced by elliptic 

equations in dimension 2. In the present context, Khasminskii's method (Stochastic Stability of Differ- 
ential Equations, Sijthoff and Noordhof,1980) leads to the study of degenerate Dirichlet problems with 
partial differential equations and nonlocal boundary conditions. 
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■ 1. Introduction 



Nonlinear oscillators subjected to vibrations represent useful models for predicting the response of me- 
chanical structures when stressed beyond the elastic limit. When the excitation is a white noise, it has 
received considerable interest over the past decades. In a previous work (1), the authors have considered 
the response of a white noise excited elasto-plastic oscillator using a stochastic variational inequality 
formulation. The results in (1) provide a framework to assess the accuracy of calculations made in the 
literature (see e.g., (4; 5; 6), and the references therein). In this paper, instead of considering white 
noise input signal whose power spectral density (PSD) is constant, we consider an excitation with a 
non-constant PSD which could be a more realistic framework. We consider the excitation as the veloc- 
ity of a "reflected" Ornstein-Uhlenbeck process. Therefore, comparing with the elasto plastic oscillator 
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excited by white noise, a third process occurs in the variational inequahty. Consider w{t) and w{t), two 
independent Wiener processes and x(t) an Ornstein-Ulhenbeck reflected process 

dx{t) = -ax{t)dt + dw{t) + l{4r)=_z.}d^/ - i{x(t)=L}d^?- (l-l) 

We have — L < x{t) < L. In this model, the stochastic excitation is given by 

-I5x{t)dt + dw{t). 

The stochastic variational inequality model is given by: 

' dx{t) = -ax{t)dt + dw{t) + l{x{t)=-L}d^^ - lw^)=i}d^^^ 

dy{t) = -{Px{t) + coy{t)+kz{t))dt + dw{t), 
< {dz{t)-y{t)dt)iC-z{t))^0, (1..2) 

When j3 7^ 0, x{t) is involved in the dynamic of y{t) and then this model will be referred as the 2d case . 
Whereas if j3 = 0, x{t) is not involved in the dynamic of ^(f) and then {y{t),z{t)) satisfy the elasto- 
plastic oscillator problem of (1) which will be referred as the Id case . 

Notation 1 .. 1 Introduce the operators 

Au:=^Uyy + ^Uxx- ccxux - (/3x + coy + kz)uy + yuz, 
1 1 

B+u := -Uyy + -Uxx- axux- {ISx+coy + kY)uy, 
B-u := ^Uyy + ]-Uxx - axux - {^x+coy - kY)uy. 



The infinitesimal generator of the process {x{t),y{t),z{t)), denoted by A is given by: 



A: (j) i-j> 




if z&]-Y,Y[, 
ifz = ±Y,±y>0. 



Notation 1..2 

0" := {-L,L) X K X ^+ := {-L,L) x (0,+oo) x {Y}; := x (-oo,0) x {-7}. 

As the main result of the paper we prove the following: 
Theorem 1..3 There exists one and only one probability measure V on &~ U 0^ satisfying 

/ A(^>dv+ / B-^dv+ / fi+(^)dv = 0, V(^> smooth. 



Degenerate DiricMet problems Related to the Ergodic property of the elasto-plastic oscillator excited by a filtered white noise3 of 46 



Moreover, v has a probability density function m such that 

/ m{x,y,z)AxAydiZ+ / m{x,y,Y)dxAy + / m{x,y,—Y)AxAy=\, 

where 

• \m(x,y,z), (jc,3',z) € ^} is the elastic component, 

• {ot(x, J,!'), (x,)') e is the positive plastic component, 

• and {m{x,y, —Y), {x,y) G is the negative plastic component. 

In addition, m satisfies in the sense of distributions the following equation in G, 

d , , d , X 1 dm I d^m 1 d^m 



and on the boundary 



ym + -[xm] + ^^m + coy + kY)m] + -^ + -^=0, m ^+ 



2m 1 ^2., 



_,^ + _[,^] + _[(/3, + eoy-fcyH + -^ + -^=0, in^- 

m = 0, in (-L,L) x (-oo,0) x {Y} U (-L,L) x (0,oo) x {-Y}. 

The proof will be based on solving a sequence of interior and exterior Dirichlet problems, which are 

interesting in themselves. We will put in parallel the Id and 2d cases, in order to facilitate the reader's 
work. In the Id case, the variable x disappears (j3 = 0), we will still use the notation A, for the 
operators defined above without x. 

Let us mention that our study presents a mathematical interest because it generalizes the method 
proposed by the first author and J. Turi (1) in the case of higher dimension. Non-local boundary con- 
ditions expressed in the form of differential equations in dimension 1 are replaced by elliptic partial 
differential equations (PDEs) in dimension 2. In the first case, there are two semi-explicit solutions. 
Thus, the non-local boundary conditions are reduced to two unknown numbers. In the second case, we 
do not know explicit formulas for solutions of the both elliptic on the boundary. In this context, these 
two solutions depend on two unknown functions respectively defined on the set (— L,L). 

In addition, the choice of the excitation (1..1) is also motivated by two technical considerations: 

• the first is to force x{[) (through the processes ^\^^) to evolve in the compact set [—L,L]. Thus, 
as part of our proof, a compactness argument allows to show the ergodic property of the triple 
{x{t),y{t),z{t ))- Note that this is not a problem in terms of applications, because if we choose L 
large enough, then the process x{t) is similar to an Ornstein-Ulhenbeck process. 

• the second is the uncorrelation of w{t) and w{t). In our approach, based on PDEs associated to 
the triple {x{t),y{t),z{t)), we avoid the appearance of cross-derivative terms in the infinitesimal 
generator A. In the case where w{t) and w{t) are correlated, these cross-derivative terms yield 
more technical difficulties. 
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2. The interior Dirichlet problem 

In this section, we prove existence and uniqueness to the homogeneous interior Dirichlet problem. 

2.1. Some background on the interior Dirichlet problem in the Id case 
Let us recall the interior Dirichlet problem from (1). Let y\ > 0, 
Notation 2..1 

Di := i-yuyi) X i-Yjy, D+ := {OJ,) X {¥}■ := {-yuO) x {-¥} 

and 

Df := {-yi + e,yi - e) x {-Y,Y),e > 0. 

Denote fi := inf{f >0,\y{t)\=yi} and consider (j) G L°°{—Y,Y). We will use the following notation 
Ep(-) := E{- |(y(0),z(0)) = p}. It is shown that E(y_^)(0(z(Ti))) solves a nonlocal Dirichlet problem: 
Find Tj e L~(£)i) n '^{DDye > such that 

ATj=OinDi, B+Tj=OinDf, B_Tj=OinZ)]" (2..1) 

with 

77(yi,z) = (/)(z), T7(-yi,z)=0, if ze{-Y,Y). 

Since T7(yi,7) = <^(y) and ri{—yi,—Y) = 0, there are semi-explicit solutions by solving ordinary dif- 
ferential equations for rj on the boundary at z = 7, and z = —Y respectively, 

i]{y,Y)=r]Yl{y,n) + (j){Y)I{0,y), 0<y^n; riiy,-Y) = r]-Yl{-yi,y), -yi^y<0, 
where Tjy and Tj_y are constants, and 

Xfexp(coA2-K2yt7A)dA 



I{a,b) := 



/o^'exp(coA2-K2m)dA 



The nonlocal condition is restricted to the value of these two constants. Based on these semi-explicit 
expressions a subset K of H^{Di) is defined for proving existence of the solution to (2.. 1) , by 

vi-yuz) = ^{z), v(->"i,z)=0, 
K=( v(y,7)=vr/(y,yi) + (/)(mO,y), O^y^yi 
v{y,~Y) = v-Yi{-yi,y), -yi^y^O 
. vy,v^y are constant with jviyl < II^Hl^- 

The set K is convex and not empty if <j){z) & (—7,7). We take v(y,z) = ^(z)/(0,>')l{j,>o}- 
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2.2. The interior Dirichlet problem in the 2d case 
Notation 2.. 2 

A, :=(-L,L)x(-yi,yi)x(-y,y), 
A+ := (-L,L) X (0,yi) x {7}, A^ := (-L,L) x (-yi,0) x {-7} 

and 

Al := {-L,L) X {-yi + e,yi - e) x (-y,y),Ve > 0. 

Denote fi := inf{r > 0, |>'(f)| =y\} and consider (j) e L°°{{-L,L) x (-7,7)). Similarly as before 
we use the notation Ep(-) := E{- |(x(0),.y(0),z(0)) = p}. We want to define E(;f^j,_^)(^(x(fi),z(Ti))) as 
the solution of the interior Dirichlet problem stated below. 

Problem 1 Find rj e L°°{Ai)n'^^{Af),ye > such that 

AT7=0in4i, S+Tj=0in4j", S_Tj=Oin/if 

and 

nxi±L,y,z) =0 in (y,^) e x (-7,7), 

T7(x,yi,z) = <?>(x,z) in (x,z) e (-L,L) x (-7,7), 
77(x,-yi,z) = in (x,z) e (-L,L) x (-7,7). 

This is formal. We should consider the case of (j) smooth first and precise the functional space, then 
proceed with the regularization. 

As in the Id-case, this problem is a nonlocal problem but the boundary condition are in two dimen- 
sions. Thus we need to solve partial differential equations for 77 on the boundary at z = 7, and z = —7 
respectively. Here we do not have semi-explicit solution, indeed tj(x,>', 7) solves 

Z?+T7=0on(-L,L)x(0,yi)withT];,(±L,y,7)=0, ri{x,yi,Y) = (j){x,z) (2.2) 

with 77(x,0,7) = r]Y{x), and r]{x,y,-Y) solves 

B_T/=0on(-L,L)x(-j;i,0)withT/;,(±Z,,);,-7)=0, T/ (x, -y 1 , -7) = (2..3) 

with TJ (x, 0, —7) = 77-y (x) where t]y (x) and 777 (x) are unknown function with || Tjiy IIl"" < || ^ ||z,~ • Next, 
we give a convenient formulation of the boundary condition (2..2)-(2..3). 

2.2.. 1 Boundary conditions 

In order to reformulate Problem 1, we consider first the equation (2.. 2) on the boundary A^. Define 
J3"'"(x,j) the solution of the mixed Dirichlet-Neuman problem 

B+/3+ = 0, in(-L,L)x(0,3;i), 

li + {±L,y) = 0, in(0,yi), 

(2..4) 

j3 + (x,yi) = 0(x,7), in (-L,L), 
j3 + (x,0) = 0, in (-L,L). 
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Proposition 2.3 If 0(x,r) eH^{-L,L) then there exists a unique solution to the equation (2..4) 

l3+eH\{-L,L)x{0,yi)) satisfying ||J3+||l=" < 
Proof. Define Dy := {—L,L) x (0,yi) and consider on //' (D^) the biUnear form 

For X sufficiently large (i^ , + ^ , Z ) is coercive. Now, define the convex set 
Kr = eH\D+),^{x,n) = <^{x,Y) and ^(x,0) = 0} 

which is not empty since e //2 {-L,L). Take z€ Ky with |jz||z,°° ^ we define to be 

the unique solution of 

bY{^l,X-^l)+H^l,X-^x)>Hz,X-^l), ^x&Kyyx&Kr (2..5) 

We have defined a map Tx {£) = ^x from Ky ^ Ky. Let us check that 

Ux\\l-^\\<I>\\l-=C. (2..6) 

Indeed, in (2..5) we take X = ^X-{^X- C)"*" e -fiTy. Hence 

-bY{^X,i^X - C)+) - A(^;t , (^X - C)+) > - C)+) 

and 

Since z + C > it follows that {^x - C)+ = 0, hence ^ C. Similarly, we check that {-^x-Cy = 0' 
hence we have (2. .6). Consider then the sequence defined by 

^7(^r\z - ^r') + - ^r^) > ^(^l;^ - (2-7) 

with 

We can take <^^(x,y) = j^^{x,Y). From (2. .6), we have ||<^_J||^„o(^q+j ^ C and from (2. .7), we have 
II^A IIhI(z)J) ^ Then, we can consider a subsequence, also denoted by ^£ such that 

^- ^ in //I (Dj) weakly and in L°°(Z)f ) weakly ★ 

also, 

^^^^ in l2(D+) strongly. 

From (2.. 7) we obtain 

feK(^,;t-^)>0, ^eKy^X^Ky. 
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We conclude easily that is a solution of (2.. 4) also ||^ ^ Il0lk°°- Now, in order to prove uniqueness, 
we must prove that a solution e //' (Dy) satisfying 

5+^=0 in (-L,L) X (0,yi) 
^,i-L,y) = UL,y)=0 (2..8) 
^(x,j;i) = ^(^,0) = 

is identically zero. 

Consider X '■= then we have 

B+x + ax + l5^y = 0mD+, 
X{-L,y)=x{L,y)=0, 
X{x,yi)=X{x,0)=0, 



hence X (^y )• This impUes 



From equation (2.. 8) we deduce ^yy e L^{Dy ). Hence ^ € H^{Dy ). In particular, ^ is continuous on 
Z)+. We have 11^11/,- = 0, so ^ =0. □ 
Now, consider the following convex sets, 

^i',M:={;t+e//i((-L,L)x(0,3;i)) satisfying (2..10); X^{x,yi)=0; \\x^\\l- ^M} (2..9) 
where 

I-L Jo ' { 5 iX^Vfx + X^Wy) + iccxx+ + {Px+coy + kY)x+)w}dxdy = 0, 
y\l/€Hm-L,L) X (O.yi)) with v^(x,0) = \i/{x,yi)=0. ^^"^^^ 

and 

K-Y,M := {X' G h\{-L,L) X (-yi,0)) satisfying (2.. 12); X'ix-yi) = 0; \\x'\\l- ^M} 

(2..11) 



where 

!\ I-y, { \ iXx ¥x + Xy Vy) + iaxx-^ + + coy - m%y ) V^ld^dy = 0, 
Vy/- e H^{{-L,L) X (-ji,0)) with v/-(x,0) = y^{x,-yi) = 0. 

These sets are not empty since they contain 0. 

Remark 2..1 • ^n{y) function of y such that n{0) = , 

X'^Kym} and {x;r; Z e ^-f,m} are //^ -bounded. 

• If X e ^^'.ll^lli-. denoting « := j3+ +X, we have 

B+O) and max(|aj(x,0)|, |co(x,yi)|) < ||<;>||l~, 
so a maximum principle implies ||ffl||L"' < 



(2..12) 
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Using the sets and /iT^y n^n, the following result gives a convenient formulation of the bound- 

ary conditions in Problem 1 . 

Proposition 2. A The Problem 1 can be reformulated in the following way: find 77 e Lr{Ai) n 
^'^(4f),Ve>0 such that 



ATj=0in4i, ri{x,y,Y)-P+{x,y)eKY^yii, riix,y,-Y) e K_Y^y\i 



and 



r],{±L,y,z) = 0, in {y,z) G {-yuyi) x {-Y,Y), 

ri{x,yi,z) = <l>{x,z), in e (-L,L) x 

'n{x,-yi,z) = 0, in (x,z) € (-L,L) x (-7,7). 

Proof. First, we can obtain the generic solution of (2. .2) by considering any function which satisfies 

X+ GH\{-L,L)x{0,n)) 

and 

B+X^ = 0, x;t{±L,y)=0, z+(x,yi)=0. (2..13) 

Note that we have not defined the value of X'^ for y = 0, hence X'^ is certainly not unique. We add the 
condition that X'^ is bounded by || <j) ||/,~. We define in a similar way the function x~ such that 

X-€H\{-L,L)x{-yuO)) 

and 

B-X'=0, x.7(±L,y)=0, X~{x-yi)=0. (2..14) 

Then, interpreting (2. .13) as (2.. 10) and (2. .14) as (2. .12) repectively, we obtain j+ e and;|f~ S 

/r_y |||^||. Hence, the set of solutions of (2. .2) and (2. .3) can be written as follows 

n{x,y,Y)-li^{x,y) eKY^ll^ll,y>Q; 77(x,>',-7) G /r_yj|^||, y < 0. 

□ 

2.2..2 Approximation (part 1) 

We study Problem 1 by a regularization method in the next proposition. Define : = A + 1 ^ . 
Proposition 2..5 The following problem: find 77*^ e L°°(zii) n//'(4i) such that 1177*^11/,- < 

A^?7^=0inAi, r]^{x,y,Y)-P+{x,y)GKY^\\^, r]^{x,y,-Y) G K_y,\\^ (2..15) 



and 



r]l{±L,y,z) = 0, in (y,z) G x (-7,7), 

77,^(x,>',±7) = 0, in (x,Ty)e(-L,L)x(0,yi), 

■q^xJuz) = <l>{x,z), in (x,z) G (-L,L) x (-7,7), 

7]'{x,-yi,z) = 0, in (x,z) G (-L,L) X (-7,7), 



has a unique solution. 
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As in the Id case, we formulate a variational inequality to prove existence of solutions. In the present 
context, we consider a convex subset of H^{Ai) which is adapted to the two dimensional boundary 
condition by 

W{x,yuz) = ^(x,z), for (x,z) G (-L,L) x (-7,7), 
K=l v/(x,-yi,z)=0, for (x,z)e(-L,L)x (-7,7), 
VA(.,.,7)-j3 + (x,3;)€/:y,[|^||^^, 

^ v^(.,.,-7)e/r_^jl^ll^„. 

Proposition 2.. 6 The set is a closed non-empty subset of H^(A\). 

Proof. The fact that K is closed follows from the continuity of the trace operator. Now, pick the 
function, 



W{x,y,z) = |- [^{x,z) - ^^{x,Y) - \^{x,Y)\ l{y^o} + (^ + ^ )^+(^,y) 1(3.^0} • 



'27 



(2.. 16) 



We have 



W{x,yi,z) 
Wix, -yi,z) 

¥{x,y,Y) 
Vix,y,-Y) 



<l>{x,z), 
0, 

j3+(x,y), >'>0, 
0, y<0. 



So, if (l){x,z) eH\dAi), (|)(x,7) e//i(-L,L), then the function defined by (2.. 14) belongs to A". □ 

Remark 2..2 IfueKandwGH^ {Ai) with w{x, ±yi,z) = and w{x,y, ±7) = 0, for < < yi then 

u + w G K. 



Consider the biUnear form 



a{u,v) = - {euzVz + UyVy + UxVx}dxdydz 

^ JAi 

+ / {Px + coy + kz)uyvdxdydz— / yuzvdxdydz+ / axuxvdxdydz. 

JAi JAi JAi 

Equation (2. .15) is formulated as follows a(M,v — m) > 0, & K,u £ K. 
Proof of Proposition 2. .5. First, existence is proved by variational argument. For X sufficiently large 
a{u,v) + X{u,v) is coercive onH^{Ai) and for / e we can solve the variational inequality 

a{u,v — u)+X{u,v — u) > (/,v — m), G K,u G K. 

We define the map m = 7\w where u is the unique solution of 

a{u,v — u) + X{u,v — u) ^ A(w,v — m), \/v GK,uGK. 

The following lemma shows that Tx is a contracting map. (see proof in Appendix) 
Lemma 2..1 If < ||0||z,~ then ||m||£~ ^ \\(i>\\L-. 
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Moreover, taking v = mq € A" we deduce 

a{u,u) + X |m|^2 ^ a{u, mq) + A (m,mo) + A7|mo — • 

That implies 



ll"llffi(4i) 



for a constant M which depends only of A , e, norm of mq and 7. 
Now, we define 

^ = {w e : \\w\\l- ^ UWl-, Ikll^i < M}. 

We have 7;^ ^ ^'^(^i) ^ L^iM) is continuous, maps ^ into itself and ^ is a compact subset of I?. Then 
Schauder's theorem imphes Tx has a fixed point u€K which satisfies 

a{u,v — u)^0, yvGK. 

Let us check that u is solution of (2.. 15). As uGK, we have 

B+u = in 4 + 

B_M = in 

u{x,yi,z) = in (-L,L) x 

M(x,-yi,z) = in (-L,L) X (-y,y) 



Moreover, 



Vv e =^ := |v e //^ (4i) such that | 



v{x,±yi,z) =0 
v(x,j,±7)=0, 0<±j<>'i 



we have u + vGK and a{u,v) > 0. Then, 

A^M = in the sense Jif' 

and integration by parts gives 



2 



/ u^ix,y,Y)v{x,y,Y)6xdy - - / / M2(x,3;,-y)v(x,j,-y)dxd); 

J-y, IJ-LJQ 

/ / Ux{L,y,z)v{L,y,z)diydiZ- / / Ux{-L,y,z)v{-L,y,z)diydiZ = 0. 
J— ft J—Y J—yt J—Y 



-yiJ-Y J-yi- 
Uniqueness of the solution to problem (2.. 15) comes from HmUl^ ^ ||<?'||l"'- D 
2.2.. 3 Approximation (part 2) 

When is smooth, we can exhibit a solution to the problem (2.. 17) by extracting a converging sub- 
sequence of Tj^. Let 6 a smooth fonction such that d{±.yi) = 0. Denote niy) := y^O{yY for some 
PA- 
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Proposition 2..7 The following problem: find 7} e L°°{Ai) such that Tja; e (4i), 

ATj=0in4i, 'n{x,y,Y)-l3+{x,y)eKY^\\^, ■n{x,y,-Y) e K_Y,im (^-l'^) 



and 

rix{±L,y,z) =0 in (y^z) € (-yiji) x (-Y,Y), 

ri{x,yi,z) = 0(x,z) in G (-L,L) x (-7,7), 

?7(x,-yi,z) = in (x,z) e (-L,L) X (-7,7), 

has a unique solution. 

As in the Id case, the key ingredient of the proof is to boimd uniformly the norm of first derivative 
w.r.t. z using the auxihary function n. 

Proof. From the previous section, we have i|?7'^||oc < ||0||oo, hence rj'^ ^ rj in L°°*. We also have 
a{ri^ ,uo — rj^) > 0, for some mq € K. So, we deduce estimates in the following lemma: (see proof in 
Appendix) 

Lemma 2..2 We have 

e / in^fdxdydz < C; / {ri^)^6xdydz < C; / (ri^fdxdydz < C. 

JAi JAi JAi 

It is hcit to test (2.. 15) with r]^y^P~^e'^. We have 



A\ 



2 < + \ nyy + \n^c+ yn! - ' - iP^ + coy + fe) n') n^y^"' 'e'^ = o. 



So, we obtain 

3 fo ff{{ri^{x,y,Y)f + {n'Ax,y,Y)WP-'d^dxdy 

+ 3 lol'n {('7'(-^'>'' + in^x^y, -Y)f}\y\'P-'d^dxdy 

+ jIa, ri^ri^{ii2p-l)y^''-^9^i+y^P-^2qe'e^'i-^)+2{Px + coy + kz)y^P-^e^i)dxdydz 
+ Ia^ nx^z {ocxy^''~ ' d^i)dxdydz. 

(2.. 18) 

Moreover, Remark 2..1 allows to bound the two integrals on the boundary that yields the following 
estimate: 

Denoting v*^ := tj^tt, we have || v*^ j ^ C;i so we can extract a weakly converging subsequence 
v*^ ^ V in H^{Ai) and v = rfn € H^{Ai). We can check that tj satisfies the boundary condition of 
Problem 1 . First, let us check that 

7i{y)Ai]=0mH-\Ai), n{y)r],{±L,y,z) =0m {m{{-L,L) x (-ji,yi)))'. 
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As e H^{A\) and kA^t]^ = 0, we have 

—A^v^ = f{r] ,r]y) in strong sense 

with/(Tj,Tjj,) := -\{n"ri^ + 2n'T]^} + {coy + kz + ax)n'ri^ . 
We obtain that V(|) e H\A^),^{x,y, ±Y) = and (^(x, ±yi,z) = 0, 

Now, when e goes to 0, we have 

t; f Vy'i'y + l; ( Vx<i)x+ j {axVx + {fix + coy + kz)Vy-yv^)<i)=[ f{r],r]y)<ii. 

Z JAi Z J Ax JA\ JAi 

We deduce we have in (zii), firstly —Av = f{ri,riy) which is equivalent to ttAtj = and secondly 

that choice of test function implies nrix{±L,y,z) =Oin {Hi{{—L,L) x (— 
Then, we check that 

■n{x,y,Y)-l5+{x,y)eKY^llpf t/(x,);,-F) e 
We know that tj^ e i/' (4] ), its trace is well defined and satisfies 

7(77^)(x,y,F)=z+'' + j3+; X^^' & Ky.Uf y>0 

with 

iiz*''IIl-s;ii</>iIl- (2..19) 

We also have X '^jX^'^ satisfy respectively (2. .20) and (2. .21) 

f-L fi' { \ + X^^'Wy) + {axy^'" + + c^y + A;y);t+'^) v^jdxdy = 0, 

\I\^^H^{{-L,L) X (0,j;i)) with V^(x,0) = v^(x,j;i) =0. 



(2..20) 



and 



' V(. + ;tj ' V^) + (a^;C^ '^ + (j3x+coy-A:y);i:5, '^)vA}dxd3; = 0, 
\/\^&H^{{-L,L) X (-yi,0)) with va(x,0) = VA(jc,-j;i) = 0. 



(2..21) 



First, we study convergence of the sequence x"^'^ and we deduce the PDEs satisfied by limg-^o/f^'^- 
In particular (2. .19) implies 

;C+'e in L^{{^L,L) x (0,yi)) weakly, 

^-,e _^ ^~ mL^{{-L,L) x {-yi,Q)) weakly. 
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And (2..20) and (2..21) imply 

< C and X'^'^^ X'^^ in weakly, 
ll/f '^^llffi ^ ^ X~'^^ ^ ;t;~;r in weakly. 

Denote 1^=^'^ ■=X^-^y^. From (2..20), we obtainB+,^+'« = in {~L,L) x {-yi,yi), in H-^{{-L,L) x 
(0,yi)). Since the operator B+ is strictly elliptic then (*+''^ e H^{{—L,L) x (0,yi)). We also have 
^;^'%±L,y) = in {H^ {Oji))'. As ^+'^ € H^{Ai) and 7rfi+x+'^ = 0, we obtain 

= g{x'^'^,Xy'^) in a strong sense 

with := - 2);(ax + coy + A;F)} - 2);;t+'^ We obtain that Vy^ e H\{-L,L) x 

(Oji)),v^(x,0) = VA(x,:yi)=0, 



Now, when e goes to 0, we have 
rL j-n 1 



/ r ^{^x'i'x+^yWy) + {ax^x+{^x+coy+m^^-y^z)'i'= I r S{X,Xy)^- 

J — LJO ^ J — L JO 

We deduce that in //~^((-L,L) x {0,y\)), we firstly have = 8{X~^,X^)> which is equivalent to 

Dndly that choice of test functions implies y'^Xxi^^^y) 

summarize, we have 
and 



y^B+x'^ = and secondly that choice of test functions implies y^Xxi^^^y) = in {{0,yi)))'. To 



Hence 
and 



-B+^+ = 8{X^,X^) in {-L,L) x{0,yi), 
^ + {±L,y) = in {0,n), 

^+{x,yi) = in (-L,L). 

nx+eH\{-L,L)x{0,y,)), U^l- ^ Ml- 



nB+x^ 



in {-L,L) X (Oji), 



Similarly, we have 
and 



Kxt{±L,y) = in (0,3;i), 
X+{x,yi) = in {-L,L). 



nx-eH\{-L,L) X (-yi,0)), Wx'Wl- < WHl- 

nB+x~ = in (-L,L) x (-ji,0), 
^Xx{±L,y) = in (-yi,0), 
^X'{x,-yi) = in (-L,L). 



(2..22) 



(2..23) 
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First, 7(?r?7'^) ^{x^ + in H\{—L,L) x (0,>'i)) weakly. Secondly, the weak convergence of 
nrj^ nrj in H^{Ai) implies the weak convergence of Y^Tirj'^) 7(?r?7) inH^ (dAi ). By uniqueness 
of the limit, we deduce 7(^7]) = 7i{x'^ + P'^)- Finally, we verify that 

ri{x,yi,z) =(l){x,z); ri{x,yi,z) =0. 

Using Green formula, we obtain 

vv^e//i(/ii), / / n'wy= I (t>yifn(y)da. 

JAi JAi JdAi 

Now, we can let e tend to 0, we obtain 

JAl JAi JdAi 

□ 

2. 2. .4 Approximation (part 3) 

Now, G L°°{dAi), we introduce a sequence of function ^ 0} C H^dAi) such that ^ ^ in 

L?{dAi). We denote 77* the solution of the Problem 2. .17 with (j)'' as boundary condition. From the 
previous section we have tj* gL°°{Ai) satisfies n{y)ri'^ gH^{Ai), 

AT7^ = 0in4i, r]''{x,y,Y)-P'''+{x,y)eKY^\m, r]''{x,y,-Y) G K_y,\^^ 

and 

ri^i±L,y,z) = in (y, z) e {-y iji) x {-Y,Y), 

n'{x,yi,z) = (l)Hx,z) in {x,z) & {-L,L) x {-Y,Y), 
n''{x-yi,z) = in {x,z)e{-L,L)x{-Y,Y). 

where j3*'+(x,j) e H^{{—L,L) x (0,yi)) solves the problem (2..4) with 0*(A:,y) as boundary condition. 
Moreover ||tj^||oo < H^IU and ||j3^||oo < ||0||oo- Let us check that the sequence J3* has a limit. 

Proposition 2.. 8 We have 

j3*'+(x,3;) ^ /3+ in l2((-L,L) x (Oji)) weakly, 
nl3'''+{x,y) -s- 7r/3+ in//\(-L,L) x (0,j;i))weaMy 

and the limit /3+ solves the problem (2..4) with (l){x,Y) as boundary condition. 
Proof. We have e H\{-L,L) x (Oji)) with V^(x,0) = xi/ixji) = 0, 

^// {Px'^Wx + Py^^ Wy)+ jj (axj3i'+ + (coj + A;F + j3x)j3,^'+)vA = 0. (2..24) 

In particular, the choice of t/A = 7r()')j3^'"'" with 7r(0) = K{yi) = gives 

yj {{Kl5!^-+f + {Kl5y''+f)+ll K(y)K'{y)l5y''+l5'-+ + jj {axP','+ + {c,y+kY + Px)P'y+yp'^'+ =Q. 
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This implies 

JJ{Kl5^'+f^C^ and jj{nP'y'+f^C^. 

We deduce that we can extract a subsequence such that we have 

J3*'+ 13+ in weakly and 7rj3*'+ 7rj3+ in weakly. 

Now, denote f 7r(y)j3'^'+. We have B+f = -j3*'+(^ - {I5x + coy + kY)n') - ji^'+n' . Also Vi/A e 
H\{~L,L) X (0,yi)) with v/(x,0) = V^(x,yi) =0, 

^// (^V^^ + W+Z/ («^}^ + (co3' + A:y + /3x))^)vA (2..25) 

When k goes to +«> in (2.. 25), we obtain 

B+/3+ = 0; J3/(±L,);)=0. 

Then, from (2. .26), 

VVA GH^nH^, jj P'^'+Yyy = -jj Py'+Vy + J (l>%n{x)da, (2..26) 
we deduce taking limit when k goes to +0° 

Vv^ e //2 n/Zgi, JJ p+Yyy = -JJl3+Yy + J (l>Vyn{x)da, 
and we obtain the Dirichlet boundary condition 

^ + {x,yi) = (!>{x,Y), j3 + (x,0) = 0. 



□ 



Theorem 2.. 9 The Problem 1 has a unique solution. 
Proof. Testing Atj * with Tj ^ , we obtain 



yj{Vye)' + yj{v^ef^yjiii')\eey + yjicoy+kz+M{v')\e'y + yja9\n'f 

-\J ax{r]^efn{x)da + ^ J yd^{'n''fn{z)da. 



So, we have 



(77^0)2 <;C; Ilin^ef^C. 



Testing (Atj*^) with r]^y^P ' 6^, we deduce the following lemma 
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Lemma 2.3 We have 

+ (P- +^Xi^(Tj*0)(Tj,%)3;^-i0«-20' (2..27) 

+ Ia, {coy + kz + px) {ri^e)ir]^n)yP- ' 
+ J^^ax{r]^e){r]^7t)yP-'e'i. 

Moreover, |j77*7r||^2 < C;i. 

Denoting := ri'^n, we have Hv^^Hz/if/ii) ^ C;^ so we can extract a weakly converging subsequence 
v'^ ^ V in (Ai) and v = riK £H^{Ai). Similarly as before, we can check that rj satisfies the boundary 
condition of Problem 1 which is summarized in the following lemma, (see proof in Appendix). 

Lemma 2..4 We have 7:{y)Ari = 0m H-^{Ai), n{y)rij,{±L,y,z) = in {m{{-L,L) x (->"i,>"i)))', 
'n{x,y,y)-P^{x,y)&KY^\\^\\, ■i]{x,y,-Y)&K_Y^\\^\\, and ■i]{x,yi,z) = ^{x,z); t/(xJi,z) = 0. 

□ 

2. 2.. 5 Local regularity in the interior Dirichlet problem 

In this section, we derive local regularity properties of the function v related to the interior problem. We 
recall that v:=nr] & (Ai) and we have 

-^^xx - ^Vyy + po{x,y,z)vy + Uxv^ - yv^ = rjpi {x,y,z) + i?yP2(y) (2..28) 

where we denote 

n" , , 

Po(x,j,z) :=/3x + coj + fe, px{x,y,z):=—-p^{,x,y,z)%, p2{y) = %{y). 

Lety <yi. 
Notation 2.. 10 

:= {(x,y,z)G4i, |y-yi|>5, b + > 5} 

/ii(5,r) {(x,y,z)e4i(5), |z-i'|>r} 

:= H\-Y,Y;H\-L,L;H\-y + b,y-m 

Recall that r\^Hx{5) means tj , r\^, r\y, tj^, r\xy,r\xz,r\7y.nxsz e (5)). 
Proposition 2.. 11 We have 

V5,r>0, Tje//i(5); ^^^{Ay{b))■, T]€H^{Ai{5,r)), (2..29) 

and 

l|T?lk(5)<Mi; \\ri\\HHMSj))<M2. (2..30) 
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Moreover, the trace of tj at y = y, denoted by h{x,z) := tJ {x,y,z) satisfies 

hGH\{-L,L) X (-F,F))n^°((-L,L) x {-¥,¥)). (2..31) 

Proof. 

The proof rehes on the following estimates (see proof in Appendix). We have 
■n^n^eL^iAi); ■qj^n^eL^iAi); T/jyTT^eL^C/ii), 

For p and q large enough, we have 

□ 

3. The exterior Dirichlet problem 

In this section, we prove existence and uniqueness to the homogeneous exterior Dirichlet problem. 

3.1. Some background on the exterior Dirichlet problem in the Id case 
Notation 3..1 

Dd := X (-F,F), Du := {y,+oo) x (-F,F), 

D :=DaUDu, D+ := {y,+<^) x {¥}, D" := x {-F} 

and 

D^:=Dn{\y\ >y + e},£>0. 

Let us recall the exterior Dirichlet problem from (1). Denote f :— inf{f > 0, \y{t)\ — y} and con- 
sider h± G L~(-F,F). It is shown that E^y^){h±{z{f))) solves a Dirichlet problem: Find ^ G L~(Z)) n 
<^°(£)'=),Ve>0 such that 

A^ = OinD, B±^ = 0,mD^ (3..1) 

with 

C{±y,z)=h±iz), -Y<z<Y. 
By solving the ordinary differential equation on the boundary, there is 

^yiy,±Y)=Kexp{coy^±2kYy). 

As a bounded solution is sought, K must be 0. Hence, ^(y,±F) = h±{±Y) and then problem 3..1 was 
recast in (1) by 

A^ = OmD, ^{y,±Y)=h±{±Y), in (3..2) 

with 

C{±y,z) = h±{z), -Y<z<Y. 
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3.2. The exterior Dirichlet problem in the 2d case 
Notation 3..2 

Ad := {-L,L) X {-^-y) x (-7,7), := (-L,L) x {y,+^) x (-7,7) 
/i:=/idU4„, 4+:=(-L,L)x();,+oo)x{7}, 4" := (-L,L) x (-00,-3;) x {-7} 

and 

/if :=/i„n{b| >); + £},£ >0. 

Due to the symmetry in the exterior Dirichlet problem, we can consider positive values of y only. De- 
note T := inf{f > 0,y{t) = y} and consider h G L"((-L,L) x (-7,7)), we define ^{x,y,z) 
as the solution of the exterior Dirichlet Problem 2: 

Problem 2 Find C e Z>"'(4«) n <r''(zif ), Ve > such that 

A^=Oin/\„, £+^=Oin/i + 



and 



C.(±L,y,z) = 0, in (y,z) G (.y,+-) x (-7,7), 

C(x,y,z) = h{x,z), in (x,z) G (-L,L) x (-7,7). 



3.2..1 Boundary condition 
Find C+ e L"(/i+) such that 



5+C+=0in4 + , (3. .3) 



and 



Define 



C+(±L,y,z) =0 in (y,z) G (y,+-) x (-7,7), 

C+(x,y,z) = /z(x,z) in (x,z) G (-L,L) x (-7,7). 



hI{A+):={u:A' 



Proposition 3.. 3 Assume that there exists //eZ/jH^'*') such that// (.x,y) =h{x,Y). 
Then there exists one and only one solution to the problem (3. .3) with 

C+G//i(4 + ), K+U-^\\hi.,Y)U^. 



Proof. First, we prove uniqueness. It is sufficient to prove that if we have 1^ = 0, ^x{^L,y) = and 

£0- 
2 • 



C(x,y) = with ^ bounded then we obtain ^ = 0. Set u{x,y) := ^(x,y)exp(— ^y^) then 



11 1 

-Uxx + ^Uyy-axux-Uy{l3x + kY) + uco{-^y^-y{l5x + kY) + -) = 0, {x,y)eA+, 

Uxi±L,y) = 0, y>y. 

u{x,y) = 0, xG(— L,L). 
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We can assume that — — >'(j3x + kY) + j <0, y^yiy sufficiently large). Let us prove that ii = 0. 
Indeed if u has a positive maximum if cannot be at y = °°. But then this contradicts maximum principle 
from (3. .3). Similarily, we cannot have a negative minimum. 

Now, we adress existence. Let 

We define a bilinear form on //j ) x ^2 ) by 

f xujV , , f (l5x + coy + kY)uyV , , 

We next define 

aj{u,v) ■.= a{u,v) + Y j^^ jy^^dxdy. 
We finally define a biUnear form on //j ('^ ^ ) x ^2 ('^ ^ ) by 

HsM := '-l^^^,dxdy + \l^j^,dxdy-2l^^^^dxdy 

+''l,-(r^2d-dy + JjPx + co{^^+coy + kY)j^^dxdy 
f , , 



/X ( \ [ CQ[y-y)dy UyV 
a,,5(«,v) =a,(«,v) - j^^ ——^^-^dxdy. 

If M e ^2 (4"'"), we can compute 

= ll^j^,dxdy+y^^j^,dxdy-2l^j^dxdy 



f xu^u ^ f ilix + coy + kY) 

A. (TW '^i^- (l+y2)2 "y^dxdy 

_£0 f U 

2 Ja+ (1 + 



1 5{y-y) 4(y-y)y 



l + 5y {l + 5yy {l + 5y){l+y^) 



dxdy 
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And see that 

2 2 1 ^ 

2j4+(l+y'^)^ 2JA+{l+y^r Ja+ [l+y^Y 



f xuxU , , f (Bx + CQy + kY) , , 

/ n I 2\2 '^^'^y+ / n , 2\2 Uyudxdy 

JA+ (1 JA+ (1 +3' ) 

2 

and the right hand does not depend on 5. Moreover, we can define a constant j depending only on the 
constants a , /3 , cq , y , A:, F but not on 5 such that 

«r,5(v,v) ^ao||v||^^ij^+j, ao>0. (3. .4) 

The constant ao depends only on a,^ ,CQ,y,k,Y . If is bounded, we consider the problem 

(iy,s{u,v — u) ^ Y j ^^_|_ 2^2 ^-'^y' ^^€^r, u€K (3. .5) 

where 

K:={veH^{A+), v{x,y)=h{x,Y)} 

which is not empty from the assumption. Then from the coercivity (3. .4) and results of the theory of 
Variational Inequalities (3. .5) has one and only one solution Uys{f) (writing u for Uysif) to simphfy 
notation). If w G //j (^^) satisfies w{x,y) = then u + w ^K, hence 



f fw 
,^s{u,w) = yj^^j^^dxdy 



and thus also 



11 y — y - + 

■^Uxx- -uyy + axux + (j3x + CO ^ _^ + Coy + kY)uy + yu = y/, {x,y)eA^, 



Ux{±L,y) = 0, y€{y,+oo), 

u{x,y) = h{x,Y), xe{—L,+L). 



Also if 
then 

Moreover, 
hence 



Mf:=m3x{\\h{.,Y)\\^,\\f\\L-} 

\\Uy5{f)\\L-^Mf. 

f f(H-u) , , 

ays{u,u) ^ ays{u,H) -y [ 7^-—^dxdy 
JA'^ (,t "ry ) 
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and 

a,s{u,H) = ay{u,H)- J^^^^^-^j^—^dxdy 

^ C||«llHi(4+)ll^llwi(4H-) + 7CM^IIM-,J')lk 
where C depends only on constants a,l5,coy,k,Y. From (3.. 4), we deduce easily that 

ll«r5(/)llHi(4+)<W)- (3-6) 

Letting 5 — >^ 0, we obtain 

Uys^Uy{f) in H2{A'^) weakly and in L°°{A^) weakly-*, Uy{f)€K 
We deduce easily from (3. .5) that Uy{f) satisfies 

/(v-«) 



y{u,v- u) > y j^^ -j^^-^dxdy (3..7) 



VvGW2'(4+), v^K, u^K, IImIIoo^M/-, \\u\\^(^^) ^Cy{f) 
where again we write u for Uy{f). The solution of (3. .7) is unique. Indeed, we first claim that 

a{u,w) = j j^^j^^^^-^dxdy, VweWi(/i+), w{x,y)=Q. 

But then if m\ are two solutions 

r ^1 ^2 /" ^^1 

ay{u^y -u^) = Y f 2\2 ^xdy, ay{u^,u^ -u^) = j f——^dxdy, 

Ja+ (l+y^V Ja+ 



hence ay{u^ — u^, — u^) = 0. However, from (3. .4) we also have 

ay{y,v)^a4v\\l,^^^y Vve//](4+) 

Therefore — u^ = 0. We next consider a sequence ^" with C,^ = ||/j(.,y)||/,«» given by the solution of 

a,(C"+',v- > r (\+))2 ' dxdy, (3..8) 

where 

C"+'€K, yv€W^iA+), veK, ||C"+i||L-^^max(||/t(.,F)||i~,||C"|k). 
Considering we have 

ay{i;\v-i;^)^yjj^^dxdy. 
Take v = ^^ — (^^ — ^'')+, which is admissible, then 
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and 



from the assumptions. Therefore (^^ — ^'')+ = which implies < and also 

Then by induction if ^ ||C"lk- ^ \\h{.,Y)\\L-. We obtain C"+' ^ C". U"^'\\l- ^ \\K.,Y)\\l-. 
Also 



Hence 



The sequence ^" ^ ^ is monotone decreasing and ^" — )• ^ in //jC'^^) weakly, ^ G K, ^ 

l|/^(-,J')lk- 
Hence in the Umit, 

a(C,v-C)^0, VveW2H.i+),ve/:, ^eK, ^\\h{.,Y)\\L- (3..9) 

and ^ is solution of equation (3.. 3). □ 

3. 2.. 2 The Cauchy problem 

The exterior Dirichlet Problem 2 is equivalent to 

AC = 0, {x,y,z)ei-L.L)xiy,oo)x{-Y,Y), 
C(xJ,z) = h{x,z), e (-L.L) X (-y,F), 

C(x,y,y) = i;+{x,y), e (-L,L) X (y 

C(±L,j,z) = 0, G (y,-) X (-y,y) 



(3..10) 



with 



and 



= 0, (x,);) e (-L,L) X 

= h{x,Y), xe{-L,L), 
^+{±L,y) = 0, ye{y,^) 

||C+|k<||/t(.,F)|k~. 
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This is a Cauchy problem, with z taking place of time and {x,y) as the space variable. We write it as 



+ ^ Qc>T + ^C«-«^S.-(/3^+coy+fe)Cy) =0 (3..11) 



and using the notation 



(3..12) 



j2^{z)u{x,y) := (^Uyy + ]^Uxx-axux-{^x+coy + kz)u^ 

we obtain 

-S + ^WC = 0, {x,y)&A+,z<Y, 

i;{x,y,Y) = ^+{x,y), {x,y) € {-L,L) x {y,^), 

U±L,y,z) = 0, {y,z)€{y,^)x{-Y,Y), 

C{x,y,z) = h{x,z), ix,z)€{-L,L)x{-Y,Y). 

So it is a Cauchy problem with mixed Dirichlet-Neuman boundary conditions. We consider the space 
i/j with a new norm 

2 f I ul + Uy f 

which is not equivalent to 

The norm in L^[A^) is defined by 



i"i'^=y..(rT?F'^'^'^ii* 

We define on T/j (^^) the bihnear continuous form 

1 [ UxVx + UyVy If Uyv{l+3y^) 



2 



^iz)iu,v) := -J^^-^^^dxdy--J^^-,^^-^6xdy 

dxdy. 



[axux + (j3x + coy + kz)uy)v ^ 



Moreover, 



y{\+y^f 



^{z){u,u) = ^l^^^^^2^dy--l^-,^^-^6xdy 

{axux+ {px + coy + kz)uy) u 



f [axux + (px + coy + kz)Uy)u 
+ y{l+y^)^ '^'^^ 
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and for |z| < Y, 

where the constants ao,fco depend only on a, ^,co,k,y,Y. Therefore also, 

< ^ {z)u,u » ao||M|lH^i(4+) - b J^^ jy-^-^dxdy. 
The problem (3.. 12) is equivalent to 

-(|^,w)+^(z)(C,w)=0, ywGHl{A+), w{x,y)=0 (3..13) 

and 

C(x,y,Y) = ^^{x,y) , ^{x,y,z) = h{x,z). 

We assume that there exists 

H{z)€H\-YJ-H2{^^)) with H{z){x,y) = h{x,z), V;c,z. (3..14) 
Writing l{x,y,z) = l^{x,y,z) - H{z){x,y), we deduce 

-(§ = {f^{z){x,y),w)-s^{z){H{zlw), ^w&H\{A+), w{x,y)=Q. 

l{x,y,Y) = t;+{x,y)-H{Y){x,y) 
l{x,y,z) = 
Under this form, we obtain one and only one solution 

Ux,y,z)&L' {-Y,Y-Mfi{A^)) , §-^{x,y,z) eL^ {-Y,Y-Mfl{A^y) 

where i/j q (4 + ) denotes the subspace of H\{A^) of function which vanish at y. We now prove that 

m\L-^\\h{.,Y)\\L- (3..15) 

We will consider 

A^ ■.= {-L<x<L, y<y<R} 

for R large. We begin with an approximation of the boundary condition ^+ with C,^ solution of (we 
delete +) 

\l^R,xx+\^R,yy-axl^R,x-{^x + coy + kY)t;R^y = 0, {x,y) & {-L,L) X {y,Y), 

CRA±L,y) = 0, y&(y,<>°), 

(3. .16) 

^R{x,y) = h{x,Y), x€{-L,L), 
Cr{x,R) = 0, xe{-L,L). 

We can assume h ^ 0, otherwise we decompose h~h^ — h^. We extend t!,R{x,y) by for y > R. The 
sequence of functions l^R{x,y) is increasing and ||CR|k°° ^ ll'j(-i^)lk°°- Let 0{y) = lif0<>'<2 and 
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if y > 1 be a smooth function. We may assume y < f • Let v G {A^).v G K and test (3.. 16) with 
"(1:^^2)2^ which vanishes a.ty = y and y = R. Setting ©^(y) = we obtain, 

2/?i4+ (l+/)2 (l+/)3j 

and thus also a{^R,v6R — ^r) = 0. Recalling that 

a{u,u) + r J^^ ■jY:^dxdy > ao||M||^^i(^+j, 

we get 

«olK«ll^i(4+) < a{CR,veR)+C\\h{.,Y)\\l^, 

from which we deduce easily, 

We then consider the limit 

^« ^ ^ monotone increasing , ^r^^ in H2{A'^) weakly 

Hence, 

«(C,v)^fl(C,0 

which implies that is the solution ^+ of (3. .3). We next consider the approximation of (3.. 10) for 
y <y <R. We write 

aCr = 0, (x,y)e4'^x(-y,y), 

CRix,y,Y) = CR{x,y), (x,y)e4^, 

C«(x,y,z) = hix,z), (x,z)e(-L,L)x(-F,F), (3..18) 

C«(x,/?,z) = 0, {x,z)ei-L,L)x{-Y,Y), 

Cfi,.(±L,y,z) = 0, (y,z)€(y,R)x{-Y,Y). 

We write (3.. 18) as a Cauchy problem 

^ + s^{z)Cr = 0, in 

CRix,y,R) = CR{x,y), {x,y)eA^, 

CRA±L,y,z) = 0, {y,z)G{y,R)>^{-Y,Y), (3..19) 

= h{x,z), {x,z)€{-L,L)x{-Y,Y), 

Cr{x,R,z) = 0, {x,z)€{-L,L)x{-Y,Y) 
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and extend by for >>>/?. If w e //j {^^)jw{x,y) = 0, we can write 

-{^,wdR) + S^{z){CR,dRW) = 0, 

^R{x,y,Y) = ^^{x,y), {x,y)eA^, (3. .20) 

^R{x,y,z) = h{x,z), {x,z)€{-L,L)x{-Y,Y). 

However, from (3.. 19), we deduce 

KrWl- < \\h{.,Y)\\L-. (3..21) 
We can pass to the limit in (3. .21) and check that ^ solution of (3. .13). This proves (3. .15). 

3.3. Local regularity in the exterior Dirichlet problem 

In this section, we derive local regularity properties of the solution C, to the exterior Dirichlet problem. 
Recall that J < y\. 

Notation 3.. 4 

Ad{8) := {{x,y,z) € Ad, y<-y-8}, /i„(5) := e 4„, y>y + 8}, 

A-{5):=Ad{5)C^A-, A+ {5) := Au{5)r\A+ , 
Ad{5,Y):={{x,y,z)&Ad{5), -Y + y<z}, 4„(5,7) := e 4„(5), z<Y-y}, 

Hd{5):=H'{-Y,Y-H\-m-H\-yi-d-n + m, 
Hu{5) := H\-Y,Y;H\-L,L;H'{yi - 5,n + 5))) 

and 

A{6):=Ad{d)UA,{d), A{6,r):=Adid,y)UA,{6,r), H{d) := Hd{d)nH,{d). 

Again, thanks to the symmetry in the exterior Dirichlet problem, we consider only positive values 

of y. 

Proposition 3.. 5 [Local regularity of the exterior Dirichlet problem] We have 

V5,r>0, Ce//„(5); C e ^ GH\Au{5,r)), 0.22) 

and 

IICIk„{5)=^A^i,||^||; !ICII//2(4„(5,r))< ^2,11011- (3..23) 
Moreover, the trace of C, aty = yi, denoted g(x,z) := ^(x,y,z) satisfies 

geH\{-L,L) X (-F,F))n<^0((-L,L) x {-Y,Y)). 
Proposition 3.. 6 [Local regularity of the exterior Dirichlet problem on the boundary 2 = 7] We have 

V5>0, t;&H\A+{d)); |KIIh2(4+(5))<M3,||^||. (3..24) 

Similar results hold for negative values of y. Proofs of Proposition 3. .5 and Proposition 3. .6 rely on 
similar estimates related to the local regularity of the interior Dirichlet problem. 
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4. The ergodic operator P 
Notation 4..1 

r± := (-L,L) X {±yi} x {-¥,¥); f ± := (-L,L) x {±y} x {-Y,Y) 

and 



ri:=rfur+; r:=r-ur+ 

4.1. Construction of the operator P 

Consider := e L"(ri). Following the same procedure of the Id case, we first solve the 

interior Dirichlet problem for tj with the boundary condition 

77 = ^+ inf+, 
7] = <p- in . 

Then 

||77|k^^max(||0+||,||0_||). 
Then, we solve the exterior Dirichlet problem for ^ with the boundary condition 

C = TJ inr+, 

C = n inf-. 
Then 

For pi e Ii, let us define 

4.2. Probabilistic interpretation of the operator P 

Let us recall from Khasminskii (7) the probabilistic interpretation of the operator P. For p & t and for 
pi e Fi , we denote by 71 ; . ) the distribution of z ( Ti ) starting from p and by fiP 1 ; • ) the distribution of 
z{f) starting from pi . We showed 

n{p) = j.^<i>{u)7i{p;du) and i;,{p\) = j_n{v)Y{pi;dv). 

Also, using Fubini theorem we can write 



where 



P0(pi) = /_ (^{u)yyi{pudu), 
Jri 

ffi{pi,du) := j_y{p\;dv)ji{v;du). 
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By construction, the operator P is the transition probability associated to the Markov chain 

where fi_o = and 

Tk+i ■.= mf{t > b(OI = yh ^hk+i :=inf{? > T*+i, \y{t)\ =yi}. 

Note that Ty^+i = to ^ and fi^^+i = ti o 9f|^^^ where dt is the shift operator. 

4.3 . Ergodic property for P 
Theorem 4.. 2 The operator P is ergodic. 

Proof. A Borel subset of Fi can be written as B := B- x {yi} UB+ x {— >'i} with B+,B- are Borel 
subsets of {-L,L) x {-Y,Y). We have 

B+ = {{x,z) : (x,yi,z) Gfi} andfi- = {(x,z) : {x,-yi,z)&B} 

also 

1b+ (x, z) = 1b (x, y 1 , z) ; 1b_ (x, z) = 1b (x, -y i , z) • 
Consider 0+ = Ib^. and 0_ = 1b_ in the interior Dirichlet problem, then 

C(x,yi,z) if3'=3'i 



C(x,-ji,z) ify = -yi. 



P(1b)(x,);,z) = 

Let p,p Gfi, define 

XpAB)-=P{lB){p)-PilB){p). 

We will prove the ergodic property of the operator P in the four following steps. 

1 . Doob's criterion from (3). 

The operateur P is ergodic if we prove the following Doob's criterion 

supXp^p{B) < 1, Vp,p G Fi and VS. 

2. Negation of Doob's criterion. 

Suppose Doob's criterion is not verified. Then, there exists two sequences pk := {xii,ykjZk), Pk '■= 
{xk,yk,Zk) infi anda sequence of Borel subsetB^^ such that 

Denote 7]*^ and the solution of the interior and exterior Dirichlet problems with = lB^^. and 
0- = lB_;t. where B+j; and are deduced from Bj. as previously. We have 

Kk:PkiBk) = C''{xk,yk,zk) - C''{xk,yk,zk)- 

Now, extracting a subsequence of pk and pk, we deduce 

Pk p* and pfc p* in Ti. 
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Also, 
and 

From (3.. 22), we know that 



p* = {x*,yuz*) or {x*,-yuz*) 
P* = {x^,yut) or {x^,-yi,t). 



then 



IK'IU(5) and <M2,i 

C'^e//(5)c'r0(4(5)), 
C'^^C* in H{5) weakly, 



(4..1) 



and 



From (3.. 23), we know that 
so, we obtain 



IC*ll//2(4±(6))^^3,l, 



(4..2) 



C'^^C in '^(4^(5)). 
Proposition 4..3 Under the hypothesis that P is not ergodic, we have ^*{p*) = 1. 

Proof. We must have Um^t-^oo C'^iPk) = 1- 

(a) If |z*| < F or z*y* = -Yyi then (4..1) impUes lim;^^^ = C*(p*)- 

(b) If z*y* = Yyi then (4..2) implies lim;t_>oo i^''{xk,±yi,±Y) = i^*{p*). Indeed, as for all k ^ 
0, is continuous, there exists k such that |C'^(/'ot) — C'^{X(yi^,±yi,±Y)\ ^ 2~*. Also, 
limA-_,^ C'^IPctJ = limA-^.^ C''{xc7,,,±yi,±Y) = C*(p*). 

□ 

3. Contradiction. 

Suppose p* = {x*,yi,z*) and set 

Ei:=rin{x = ±L}; E2 := nn{z = Y}. 

Maximum principle for paraboUc operator appUed to ^* impUes p* € EiU E^. Then p* cannot be 
in Si because of the Neuman condition and p* cannot be in E2 because of the boundary condition 
z = Y. A similar argument yields a contradiction with p* = (x*, —yi,z*). 

4. Conclusion. 

From ergodic theory, there exists a unique probability mesure y* = (7^)7*^) on Ii and K,p >0 
such that 



Vn>0, 



</5:||0||exp(-pn). (4..3) 

□ 
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5. The nonhomogeneous interior Dirichlet problem 

Consider / a bounded function on Ai, we want to define ^{x,y,z) ( fo^ f{x{^)jy{^)iZ{s))ds) as the solution 
of tlie non-liomogeneous interior Diriciilet Problem 3: 

Problem 3 Find x G L'^{Ai)r\'rf^{Af)ye > such that 

AX + f = in /ii, 
B±X + f = in Af, 
X{x,±yi,z) = in {-L,L)x{-Y,Y), 
Xx{±L,y,z) = in {-y^j,) x {-Y,Y). 

Consider 

0{x,y,z) := expA(coA:z^ + coy^); A ^ 1. 
Theorem 5.. 1 The Problem 3 has a unique solution. Moreover, we have 

\\X\\L-^exp?i{cokY^ + coyl), 

where A depends on /. 

Proof. The uniqueness of a solution of Problem 3 is argued as the uniqueness for the homogeneous 
Dirichlet interior problem. The existence can be proven by the regularization technique used previously. 
Now, we give a L°° bound for X- We have 

0x = 0; = 

(^>j, = 2ycoA^; <^j,y = 2coA(^ + (2coA)')-^(^> 
(j)^ = 2zcokX<p 

and 

y^z - ax<l)x - <^>j,(j3x + coy + kz) + ^^^y + ^<fec = -2/3coAxy(/) - 2coA/0 + cqA^ + 2(coAy)^ 

= coA(;)-coA2j3xy(;) + 2(co3')^A (A -!)(;» 

> coA(l-2j3Lj;i)0 + 2(coy)^A(A-l)^ 

> 11/11 

where yi can be chosen as 2yi < and A > max(l, co(i-2|?Lyi) )- 
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Using that <j){x,y,±Y) = exp A {cokY-^ ) exp A (cqJ^ ) we obtain for z = Y,0 <y <yi. 

-ax<l)x-(j)y{lix + coy + kY) + ^<l)yy + ^<j)^ = -(j)y{x,y,Y){coy + kY + l5x) + ^(l>yy{x,y,Y) 

= kco^{x,y, F) ( 1 - 2{kY + ^x)y + 2co(A - 1)/) 
> Aco0(x,3;,y)(l -2(A:y + /3L)); + 2co(A - l)y^) 

I 

2co(A-l) 

^ llfIL 



, {kY + l3Lf 



if A > max(2W J + M±£i;L). Similar estimates hold for z = -Y, -yi<y< 0. Consider m := <^ - 
then u satisfies the following inequalities 

AM^Oin/ii; B±u ^0 in (5..1) 

and the following boundary conditions 

u{x, ±yi,z) = exp(A (cofe^ + coJi) ) in (-L,L) x (-F, Y), 

u,{±L,y,z) = in(-y-i,);i)x(-y,y). ^^"^^ 

The maximum of m cannot be attained in 4 1 or on the boundaries z — ±Y,x = ±L. It can only be attained 
for y = yi ory=yi. Hence, 

u{x,y,z) ^ exp(A(coA;F^ +co>'i)) 

which implies 

X{x,y,z) > -cxp{X{cokY^+cofi)). 
Now, consider V = — (0 then v satisfies inequahties (5.3) 

Av^Oin^i; fi±v^OinA± (5..3) 

and the boundary condition (5.. 4) 

v{x,±yi,z) = -exp(A(cofe^ + co3'?)) in (-L,Z,) x (-7,7), 

v,{±L,y,z) = in(-j;i,j;i)x(-y,y). ^^""^^ 

Again, the minimum of v cannot be attained in Ai or on the boundaries z = ±Y,x = ±L. Hence 

v{x,y,z) > -exp{X{cokY^ + coyj)) 

which impUes 

X{x,y,z) < exp(A(coA;F^ + co3'i)). 
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Now, let us derive some estimates on derivatives. Let M := exp(A(cofcF^ + co/i^)). Using the test 
function mo{x,y,z) = exp(— co(3'^ +fe^))exp(— ax^) we obtain now a priori estimates on tlie partial 
derivatives of X- We test (3) with moX 

/ Xymo+ xlm = / yX^m- yx^mo + 2 fxm- l coPxx^mo 

JAi JAi JdAi,z=Y JdAi,Z=-Y JAi JAi 

< / yX^mo- / yx^mo + 2 / fx 

JdAi,z=Y,y>0 JdAi,z=-Y,y<0 JAi 

^ C(M). 

/ Xymo^CiM) ■ [ xlm^C{M). 

JAi JAi 



Hence, 

We have then. 



/ iy'Xzfmo^CiM) and / {y'Xyyfmo^CiM). 

JAi J Ai 

The function x is smooth outside a neighborhood of y = 0. □ 

6. The nonhomogeneous exterior Dirichlet problem 

Consider the function 

:= \og(y)+K, K such that login) +K> 

and the space of functions 

X^:={uY, u€L°°{Au)}. 

We want to define E(^^y^^'j(^fQ f{x{s),y{s),z{s))ds) as the solution of the nonhomogeneous exterior 
Dirichlet Problem 4: 

Problem 4 Find ^ eX~n'^0(4f ),Ve > such as 

A^+f = in Au, 

B+^+f = in A + , 

^ixj,z) = in (-L,L) X (-y,F), 

U±L,y,z) = in (y, +-) x (-7,7). 

Theorem 6. .1 The Problem 4 has a unique solution. Moreover, we have — ^^{x,y,z) ^ W{y)- 
Proof. We justify uniqueness of the solution taking / = and setting ^ 

B,=w\^, a>l. 

We show w = 0. In particular, w satisfies 

^ = yvxW'^\ ^y = wy\ir'^ + wa\ify\if'^-'^\ = Wzw" 

^XX = WxxW"' ^yy=Wyy\l/"+2Wya\l/y\l/"-'^+Wa{\l/yy\l/"-^ + {\l/yf\l/"-^{a-l)) 
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then 



- {coy + kz + I5x){wy\i/" + wa\i/y\i/"~^} - axwx^j/" + ywz v'" = 0. 
Collecting terms related to i/a", we obtain 

+{-coy-kz- Px){wy+wa^}- axwx + yw^ = 0, 
which implies for z € (— F, Y) and y > y, 

rW„ + -Wj5, + Wj,{a-!^-Co-fe-j3x} + {7;Wyy + - 



- (coj + fe + J3x) a v^j,} - axwx + yw^ = 



and for z = 7; j > y, 



- (coy + A:F + I5x)a\ify} - axwx+ywz = 



and 



w{x,y,z) = o. 

We can find a > 1 with a — 1 small so that 

Since w as y oo, a positive maximum can be attained only for a finite y*. But this is impossible 
from the equation. So w is negative.w < 0. But, w is also positive, w > 0. Hence ^ = and uniqueness 
follows. 

Now, existence is demonstrated by the following approximation 

A^^ + f = 0, in Ar, 
B+^« + / = 0,in/i+, 
^{x,y,z)=0, ^ix,R,z)=0, 
^{±L,y,z)=0. 



Using m'' = ^'^ - v^, we obtain 



Am* > 0, in Ar, 
B+M* >0, in 4^, 

u'^{x,y,z) = -vf(y), 

4{±L,y,z)=0. 
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Necessarily, ^ 0. The sequence is, monotone increasing and converges towards a solution. Let us 
show estimates on u. Suppose / ^ 0, we will show that ^ ^ ^ i/- Then, 



Define y such that ^ ^ ||/|| . We then have, setting u = ^ —y, 

Au > 0, in /i 
B+M^O, in 4+ 
u{x,y,z) = -V^(y) 

It follows that M < 0. 

We can show similarly that ^ + ^ 0. Hence, we have 

□ 

7. The operator T and the invariant measure v 

7.1. Construction of the operator T 

Consider f & L°°{0'), following the same procedure of the Id case, we first solve the interior non- 
homogeneous Dirichlet problem for x with / at the right hand side, then we solve the exterior non- 
homogeneous Dirichlet problem for ^ with / at the right hand side and ;t as a boundary condition. 
Also, for any pi S /i, we define the operator 

T/(/5i):=^(pi)- 
This defines a hnear operator from L°° (fl) in L~(Ii). 

7.2. Probabilistic interpretation of the operator T 

For any A a Borel subset of ^, consider the two following measures of occupation of A by the process 

{x{t),y{t),z{t)) starting at p € ^, namely 

v^{p;A) := Ep{J^^lA{x{s),y{s),z{s))ds) 
v^{p;A) := Ep{J^lA{x{s),y{s),z{s))ds). 



= r[^{-coy-kz-l5x)-^] 

^ r[^(-coj-fe)-^] 

< 7[^i-coy-kz)-^] 

Co 2kY 
2 Co 



Degenerate DiricMet problems Related to the Ergodic property of the elasto-plastic oscillator excited by a filtered white noise35 of 46 



We have shown 

Xip)= [ f{q)dVx{p;dq) and ^{P)= j f{q)Av^(p;Aq). 
For any p\ G r[, we have 

T(/)(^i) ■.= Ep^{ r f{x{s),y{s),z{s))ds)+ [_E,{ fix{s),y{s),z{s))ds)dYipv,dq) 
T(/) integrates / over a time cycle of the Markov chain stopped on rf, that is {x{fi^i!;),y{fi^k),z{ti^k)). 



7.3. Construction of the invariant measure V 
Now, define 

. . Jf^/5yr/(r,yi,.)y+(dr,d.)+/^^£,r/(r,-yi,.)y-(dr,d.) 
/-\/-^2"l(r,yi,5)7+(dr,d5) +/f J_^5,ri(r,-j-i,5)7,-(dr,d5) ' 

The denominator is well defined and is stricly positive. Now, we want to solve the complete problem. 

Problem 5 Find u such that UY~^ is bounded for \y\ > \y\ and 



Au + f = 


0, 


in 




B+u + f = 


0, 


in 




B-U + f = 


0, 


in 


e- 


Ux{±L,y,z} = 


0, 


in 


{x = 



Considering Problems 1, 2, 3 and 4, in the following result we use a functional analysis method to prove 
that V is the unique invariant distribution associated to the solution of the stochastic 

variational inequality (1..2). 

Theorem 7..1 The Problem 5 has a unique solution if and only if v(/) = 0. 

Proof. Similar arguments as the Id case are considered. Uniqueness is guaranteed by the ergodic 
property of the operator P. Suppose v (/) = 0. We define % the solution of the Interior non-homogeneous 
Dirichlet Problem 3 and ^ the solution of the Exterior non-homogeneous Dirichlet Problem 4. We set 

:= X and i^^ := i^, and for ^ ^ 0, we define x^^^ by 





= 0, 


in 


4i, 




= 0, 


in 






= 0, 


in 




X^^\±L,y,z) 


= 0, 


in 





with 

;t;*+i(xji,z) = i,^{x,yx,z)\ X^^\y^-y\-,i) = ^''{x,-yi,z). 
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Afterwards, we define i^/t+i by 

A^'^+i = 0, in 

B+^*+i = 0, in A+, 

^!^+\±L,y,z) = 0, in {y,+<^) x {-Y,Y) 

with 

Similarly, we define ^^"""^ by 

A^*+i = 0, in Ad, 

= 0, in 47, 

^^i(±L,y,z) = 0, in -y) x (-7,7) 



with 

That means 
Next, we set 
then 



(7..1) 



I'if. =T/ + P(T/) + ...+P*(T/). 
Now, let us remark that x'' satisfies the following equation: 

Af + f = 0, in Ai, 

B+f + f = 0, in A + , 

B-f + f = 0, in A^, 

fA±L,y,z) = 0, in {-00 -y) x {-¥,¥) 

with 

X\x,yi,z) = l^-^xji,?); x\x,-yi,z) = ^''-\x,-yi,z) 
satisfies the following equations: 

a|^ + / = 0, in Au, 
B+l'^ + f = 0, in A+, (7..2) 

^K±L,y,z) = 0, in {y,+<^) x {-Y,Y) 

^Hx,y,z)=f{x,y,z) 



Degenerate DiricMet problems Related to the Ergodic property of the elasto-plastic oscillator excited by a filtered white noise37 of 46 
and 

Al^ + f = 0, in Ad, 
B+l' + f = 0, in AJ, (7..3) 

^K±L,y,z) = 0, in {-^-y)x{-Y,Y) 

The condition v(/) = means 

/ Tf{x)dn{x) = [_ Tf{r-yus)dni{r,s)+ [_ Tf{r,yi,s)dn2{r,s)=0. 
From the estimate (4.. 3), we obtained 

I converges in L°° {fl ) . 

Now, we notice that %^ — % is a solution of the interior homogeneous Dirichlet problem with — 
%)\rx — \t\ I* — § is a solution of the Exterior homogeneous Dirichlet problem with — 
^)lf = {X^-X)\t- Then, we obtain 

We can extract a subsequence such that 

Moreover, % satisfies equation: 

= 0, in Ax, 

B+X + f = 0, in A + , 

B-X+f = 0, in A-, 

U±L,y,z) = 0, in {-y,j,) x {-Y,Y) 



(7..4) 



X{x,yi,z) = ^{x,yi,zy, X{x,-yi,z) = Ux,-y\,z) 

and I satisfies equations: 

a|+/ = 0, in Au, 
B+l+f = 0, in A+, (7..5) 

U±L,y,z)+f = 0, in {y,+oo) x {-¥,¥) 

^{x,y,z)=x{xj,z) 
Al+f = 0, in Aa, 
B-l+f = 0, in A-f, (7.6) 
U±L,y,z) = 0, in {-^-y)x{-Y,Y) 
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Ux,-y,z) = x{x,-y,z). 

Then, we must have 

X = | in >^<±>'<>'i 

Thus, the function 

jx in/ii, 
u= < ~ 

in /if 

is the solution of Problem 5. 

Now, suppose that Problem 5 has a unique solution u. Considering the same sequences and we 
have that m — is a solution of the interior homogeneous Dirichlet problem with (m — x'^) | = {u — 
|'^"')|fj andw — I ^ is a solution of the exterior homogeneous Dirichlet problem (m — I*) I j=. = {u — x'^)\f'. 
Hence, 



and so. 
We have 



lll^llL'«'(ri) ^ ^- 
I' = {k + l) / T/(x)d;r(x) + £p^(T(/- v(/))) 

n ;—f\ 



k 

Tf{x)dK{x) + 
and since the sum is bounded, we obtain 



{k+l) I T/(x)d7r(x) is bounded. 



That leads to 

/ T/(x)d7r(x) = 

That implies v(/) =0. □ 
Consider now (p a smooth function in [~L,L] x R x [— with compact support. If we take 

f{x,y,Y) = -B+cp, 
f{x,y,-Y) = -B.(p, 

then (p is solution of (5) for this /. From Theorem 7..1, we have 

I^L /-oo I^Y { ^ ^xv + ^ <P« - «^<P-i ~il^x + CQy + kz}(Py+y(p,}dv{x,y, z) 



' 2 "^'^ 2 

+ I-L lo { 2 ^yy + 2 ~ "^^"^ - {lix + coy + kY)(py}dv{x,y,Y) 



+ J!:l I-.o{^(pyy + \9xx- ax(px - {Px + coy - kY)(py}dv{x,y, -Y) = 0. 
If V has a density m, then we deduce that 

a^[xm] + -^[{Px + coy + kz)m]-y^ + ^^ + ^^=0 in (0,L) x R x (-7,7) (7..7) 
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in the sense of distributions. If we test (7. .7) with (p and integrate by parts, we obtain 



j ym{x,y,Y)(p{x,y,Y)dxdy + j ym{x,y,-Y)(p{x,y,-Y)dxdy 

iL lY J J 



and comparing with (7. .7) 

rL rO pL r+oo 

/ / ym{x,y,Y)(p{x,y,Y)dxdy + / / ym{x,y,-Y)(p{x,y,-Y)dxdy 
J-lJ-oc Jo Jo 

i-L roa 2 J 

- J ^J^ m{x,y,Y){-(pyy + -(p„-ax(px-{Px+coy + kY)(py+y(p}6xdy 
J m{x,y,-Y){-(pyy + -(p^-ax(px-{l5x + coy-kY)(py-y(p}dxdy = 



, r 1 1 

we finally deduce 



ym+-^[xm] + -^[{Px + coy + kY)m] + -^ + --^=0, on ^+ 

d r , d ,„N 1 1 d^m 1 d^m „ ^ 

-ym+-[xm] + ^^[{^x + coy-kY)m] + -^ + -^=^, on 

OT = 0, on (-L,L) X (-°o,0) X {7} U (-L,L) x (0,°o) x {-7} 
The proof of the main result is complete. 

8. Appendix: Proofs of lemmas 2..1,2..2,2..3,2..4 

Proof of Lemma 2..1. Denote 7= Notice we have 

(M-7)+(x,±yi,z) = ifxe(-L,L), zG(-F,F) 

{u-y)+{x,y,±Y) = ifxe(-L,L), 0<±);<j;i 

and 

(M + r)-(x,±yi,z) = ifxe (-L,L), ze{-Y,Y) 

{u + Y)~{x,y,±Y) = ifxe(-L,L), 0<±);<j;i. 

Choose V = M — (m — y)"*" e K and we obtain 

-a(M, (m - 7)+) - A (m, (m - 7)^) ^ -A (w, (m - 7)+) 
a(M,(M-7)+) + A(M,(M-7)+) < A(w,(m-7)+). 
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Switching the first argument by (m — 7) in the previous inequaUty, we obtain 

a{{u - r)+, {u - 7)+) + A |(m - 7)+|^2 < A (w - 7, (m - 7)+) 

and w ^ 7 imphes {u — 7)+ = 0. Taking v = m + (m + 7)" similar arguments yields (m + 7)" = 0. □ 
Proof of Lemma 2. .2. We have the following expressions, 

and 

f e 1 1 

a(Tj^Mo) = ^ri!uoz + ^riyUoy+ -riluox + {Px+kz + coy)uori^ -yuon! + «xuonx- 

Inequahty a(?7'^,T7'^) ^ a{ri^,uo) means 



{Px + coy + kz)ri^ri^ + ^j^ (j3x + coy + fe)Moi]y 
+1 [ VxUOx+ [ ax-q^-q^- [ axuor]^ 



lAi JAi 

We apply Cauchy-Schwartz inequaUty to the first two terms, then we apply Green formula to the last 
one and we get 



r r blJxJy+II^IK/ bMozl)+ t r \y\Wo{x,y,Y)\ 

J—L J —V] J A} J—L J — vi 



\y\\m{x,y,-Y)\ 

LJ-yi 

with ci,C2,C3 and C4 are positive constant, we get 

\\Vlnl\\h + H\\l2 + Wx\\l2^^^c,\^^^^ 

□ 
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Proof of Lemma 2.3. We have 

J At J At J A\ J At 



lAi JAi JAi JAi 

1 

2 JAi " J Ax 



2q [ (77y^0)(T/y^-l029-2)0' 

JAi 



. . {r^^)\x,y,Y)y'P-'9''i + - / {r^^fix,y-Y)/P-'d''' 

^J-LJ-yi ^ J-L J-yi 

-{2p-i)[ in^e){n',n)yp-^e^-'-2q [ {n^e){n^n)yp-'e^-^e', 

JAx JAi 

JAi ^ JAi 

j^J—LJ—yi ^J—LJ—yi 



and 



/ icoy + kz + Mn^n^y^P-'d^" = [ {coy + kz + ^x){r]^e){r]'^n)yP-'e"- 
JAi JAi 

I axr]lnly^P-'d^ = [ ax{n',d){n',n)yP-'d^. 

J At J At 



Testing [Arf) with T]^y^P-^d^, we deduce 



J in^^f = ^JJ'\nyn^,y,Y)y'p-'e''i--JJ''_^^^^^ 

+ \ j\f]y^)\^^y^'^)y'''~'^^~\ l[f\ri',f{x,y-Y)y'P-'d'^ 

+ / (coj + fe + /3x)(T/*0)(T/,^;r)/-i0«-i 

JAi 

+ I axir]^e){r]^n)yP-'0i. 

J At 



Proof of Lemma 2. .4. 
1. With 



:= - {cQy + kz + ax)n') - ri^n', 



□ 
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wehave'ixj/ eH^{Ai),\i/{x,y,±Y) =0 and\i/{x,±yi,z) = 0, 













/ 


/4i 


Jai 



Now, when k goes to we get 

If vy\i/y + ^ [ v;trx+ / (axv;, + (j3x+co)' + fe)v5,-);vz)v/-= / f{ri,riy)^. 

^ J A\ ^ J A\ J A\ J A\ 

We deduce that in //^'(zii) we firstly have —Av = /(i],i7y) which is equivalent to %Ar\ — and 
secondly that choice of test function implies %'r]xi^^^^y^z) = in {H^{{—L,L) x (— 

2. (a) We know that nrj^ G (^i)^ its trace is well defined and satisfies, 

Yi7tri''){x,y,Y)=K{x+''' + li'''^y, y>0 
Y{7iri''){x,y,-Y) = KX~'''; y<0 

with 

\\X^'%- < UWl- (8..1) 
and X~'''jX~^''' satisfy respectively (8. .2) and (8. .3) 

I-Llt{kiXx''Vx + X^*¥y) + {axx;t'' + {lix + coy + kY)x^'')Y}dxdy = 0, 
y\l/&H\{-L,L) X (0,yi)) with \lf{x,0) = \l/{x,yi)=0 

and 

I-L I-n { \ iXx'^Vx + XP'Vy) + {axxP" + {px + coy - kY)xP%}dxdy = 0, 
y\l/eH\{-L,L) X (-yi,0)) with \l/{x,0) = \i/{x,-yi) = 0. 



(8..2) 



(8..3) 



First, we study convergence of the sequence x^''^ and we deduce PDEs satisfied by lim/c-^o X^''^- 
In particular (8..1) implies 

X+'''^X'^ mL^{{-L,L) X (0,yi)) weakly, 
X~'^ ^ X~ '-^L^{{-L,L) ^ (-3'!, 0)) weakly. 

and equaUties (8. .2), (8. .3) imply 

|lX+'*^7r||^i «;C and ;t:+'*7r -l> ;t:+;r in H'((-L,L) x (0,>"i)) weakly, 
<C and ^ ;t":;r in //'((-L,L) x (-yi,0)) weakly. 

Denote^*'* := ;i:±V and := -X^''{^-M(^x+CQy + kY)] -lyXy''' . From 

(8..2), we have fiy^+'^ = g{x+^\Xy'^) in {-L,L) x (->^i,>^i), in H-\{-L,L) x (0,j;i)). 
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The operator B+ is strictly elliptic then (^+'* e H^{{-L,L) x (0, >"i )). We have {±L,y) = 
Oin (i/5(0,yi))'. As ^+'* e H^{Ai) and ;rB+x+'* = 0, we have 

-8+^+''' = g{x'^''',Xy''') in a strong sense 

We obtain that Vv/^eff^((-L,L) x (Oji)), va(x,0) = \ir{x,yi) =<d, 

-lJo 



Now, when k goes to 0, we have 

J-L Jo' \ "^^ ^ "^"^ ^ ^"'"''^^ ^ ^'^'^ ^ ^ ^^^''^^ ~ \-L r ^ ' 

We deduce that in ((— L,L) x (Oj^i)) we have — = which is equivalent 

to >'^B+X+ and that choice of test function implies >'2x/(±L,>') = in (//2 ((0,yi)))'. 
(b) Firstly, y{n'q'') -> 7r(x^ + j3+) in H^{{-L,L) x (0,yi)) weakly. Secondly, weak conver- 
gence of ;rTj^ — )• ;rTj in //^ {Ai ) implies weak convergence of 7(;rTj^) 7(^:77 ) in //2 {dAi ) . 
By uniqueness of the limit, we have Yi^'H) = ^iX'^ + J^"*")- 

3. Using Green formula 

Now, when k goes to °°, we obtain 

Vv/-e//o(4i)nj/^(4i), ^ ^^yy^ ^^^^^ La ^^y*^^- 



□ 



Lemma 8.. 1 We have 

Proof. Denote w:= Vx- Deriving equation (2. .28) with respect to x, we obtain equality (8. .4) 
- ^Wxx - ^Wyy +A{x,y,z)Wy + Uxwx - yw^ = -jivy - avx + pirix + (pi)^:'? + p2{y)'nxy 
Testing equality (8. .4) with n^, we obtain 



axwxwn 

4i 



- / {wx7t)^ + - I (wyTT)^ = — / Wyitn'w— [ A{x,y,z)Wywn^ — [ 
2 2 Jai Jai Jai J a 

+ / ywzWK^ — 15 / Vywn^ — a / VxW7t^+ / piTj. 

JAi J A\ JAi J Ai 

+ / {P\)xriwn'^+ I rixyP2(y)wn^ 

JAi JAi 
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which means 



- / {wx^)^ + [ {wyTty^ = — I WyJtn'w— [ A{x,y,z)wywn^ — f axwxwn 

2 J Ai 2 J Ai JAi JAi JAi 

+ / yK'^w^n{z)da - j3 / VyWTi'^ -a v 

JdAi JAi JAi 

+ / piw^n+ / {pi)xriwn^+ / rixyp2{y)wn 

J Ai JAi JAi 

+ I {wyn)p2{y)w- / Vxn'p2{y)w. 
Ja< JA] 



It is easy to verify that Vxx^,VxyK €:I?{Ai). 

Denote w := Vy. Deriving equation (2. .28), with respect to x, we obtain equaUty (8. .5), 

- ^Wxc - ^Wjy + CQW+A{x,y,z)wy + axwx - yw^ = +P1T/3, + {pi)yr] + r]yyp2{y) + ^yp2{y)' ■ (8.. 5) 
Then, we test equahty (8. .5) with tt^w, we obtain 



/ {wx'^)^ ^ — / {wyTi)^ = — [ WyTtn'w— [ C(j{w7t)^ — [ A{x,y,z)wyW7t^ — [ axwxwn 
jAi 2 Jai Jai Jai Jai Jai 

+ / v^wn^+ / ywzwn^+ / r\ypiwn^+ / r\{pi)ywn^ 

J A\ J Ai JAi J Ai 

riyyPiwTt^ + %(p2)'?r^w 



JAi JAi 

which means 



- / iwxn)^ + - [ (wyK)^ = — ( WyKn'w— ( CQ{wn)^ — f A{x,y,z)wywn'^ 

2 J Ai 2 J Ai J Ai JAi JAi 

- I axwxW7t^+ / v^wK^+ / y{wK)^n(z)do + / rjyPiwK^ 

Ja^ JAi JdAi JAi 

+ / '?(Pi)yw;r^+ / {yvy-2r]yn' + r]n")wn+ / riy{p2)'n^w. 
Jaj Jai Jai 



Again, it is easy to verify that Vyx7t,Vyyn G L?'{A\). □ 
Lemma 8.. 2 We have 

rjxzy^Ti'^ e L^iAi); r]y,/7i^ e L^iAi). 

Proof. 
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We test equality (8. .4) with y^P^^n^w^, we obtain 



= \jj wl{x,y,Y)y^P-'K^dxdy- \jj wl{x,y-Y)y^P-'K^AxAy 

-IJ^Wyyiy'P-'K^'^w,) = y^WyW^y^P-'n^'' + J^WyW,{{2p-l)y^P-V'i + y^P-'2qn^''-'n') 
= ^ JJ w]{x,y,Y)y^P-'n^dxdy- ^ JJ wl{x,y,-Y)y^P-'n^^dxdy 
+ {wyK){w,yPn^){yP-^K'i-')+qj^ {wyK){w,yPn'i){yP-'K'i-'K'). 



We have 



^ (/JtV)' {wl(x,y,Y)+wl{x,y,Y)}y^P-'n^'^6^Ay 

- \ JJ{wlix,y-Y)+wlix,y-Y)}y^P-'K'^dxdy 

+ [ A{x,y,z){wy7t){w^P7t'^)yP-'^7t''-^ + [ ax{wj,7t){w^P7i'')yP-'^7i'^-'^ 

JAi JAi 

+ [ ^Vy{yPniw,)yP-'n^+ ( av,{yPn''w,)yP-'n^ 

JAl JAi 

-[ pMyPn^w,)yP-'n^-' + f (pi).T/(/;r«w,));f-i;r« 

JA^ JAl 

- / P2{vxyn-v,n'){yPn''w,)yP-'n''-^. 



For p > 2 and some calculations give Vxzy^Ti^ GL^{Ai). 
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Now, we test equality (8. .5) with y^P ^n^w^, we obtain 



+ {p~h[ iwyn){w,yPn'>){yP-^ni-') + q [ {wyn){w,yPn'>){yP-'n''-^n') 
+ [ cow{yPn''w,)yP-^n'' 



Jai 



+ [ A{x,y,z){wyn){w,yPn^)yP-'K^-'+ [ ax{w,K){w,yPn^)yP-'n^-' 




I v,{yPK^w,)yP-'K'^- f piivy-i^K'){yPK^w,)yP-'K'^-' 




- I (pi),T/(/;r«w,)/-i;r«- / (pi).Tj(/;r«w,)/"'^'- 

J Ai J Ai 



For p^2 and q^2, some calculations gives Vy^y^Tt^ S L^{Ai). 
Lemma 8.3 For p and q large enough, we have 



□ 



Vxyzy"^ e L2(/ii); r^yy.yPn" e L2(/ii); p{z)n^,yPK'i e L^{Ai); p{z)n^yPn'i e L^{Ai). 
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